Electromagnetically induced transparency and wide-band wavelength conversion in 
silicon nitride microdisk optomechanical resonators 
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We demonstrate optomechanically-mediated electromagnetically-induced transparency and wave- 
length conversion in silicon nitride (SisN4) microdisk resonators. Fabricated devices support whis- 
pering gallery optical modes with a quality factor (Q) of 10 6 , and radial breathing mechanical modes 
with a Q=10 4 and a resonance frequency of 625 MHz, so that the system is in the resolved sideband 
regime. Placing a strong optical control field on the red (blue) detuned sideband of the optical 
mode produces coherent interference with a resonant probe beam, inducing a transparency (absorp- 
tion) window for the probe. This is observed for multiple optical modes of the device, all of which 
couple to the same mechanical mode, and which can be widely separated in wavelength due to the 
large bandgap of SisN4. These properties are exploited to demonstrate frequency upconversion and 
downconversion of optical signals between the 1300 nm and 980 nm bands. 
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Recent demonstrations of strong radiation pressure in- 
teractions in cavity optomechanics have focused on the 
resolved sideband regime, where mechanical sidebands 
of the optical mode lie outside its linewidth jlj. Such 
systems have been used in laser cooling of a mechanical 
oscillator to its ground state @, H[ , coherent interference 
effects such as electromagnetically-induced transparency 
and absorption 0,0], parametrically-driven normal mode 
splitting [6j-|8j], and observing energy exchange between 
the optical and mechanical systems [H, [§] . In the con- 
text of chip-based micro- and nanophotonic devices, silica 
microtoroid cavities @, Q and silicon photonic/phononic 
crystal resonators (optomechanical crystals) [H, [a, El are 
two successful systems that exhibit significantly different 
parameter regimes. Silica microtoroids typically support 
ultra-high quality factor optical modes (Q > 10 7 , decay 
rate k/2tt ps 10 MHz) that are coupled to > 50 MHz fre- 
quency mechanical modes with a zero-point optomechan- 
ical coupling rate go/lir > 1 kHz. In contrast, silicon op- 
tomechanical crystals have Q « 10 6 (k/2it > 200 MHz) 
modes that couple to few GHz mechanical modes with 
.g /27r fa 1 MHz. 

Here, we consider a chip-scale platform formed by sili- 
con nitride (SiaN4) microdisks, in which Q pa 10 6 modes 
couple to a 625 MHz mechanical radial breathing mode 
with 5o/27r 8 kHz. This sideband-resolved device com- 
bines some advantageous features of the aforementioned 
systems, such as high mechanical frequency, a straightfor- 
ward optical mode structure with multiple high-Q modes 
coupled to the same mechanical mode, and operation 
in a material whose broad transparency and low nonlin- 
ear loss enables large intracavity photon numbers to be 
achieved over a range of wavelengths. We demonstrate 
optomechanically-mediated electromagnetically induced 
transparency (EIT) and wavelength conversion Ill4l4j. 



up- and downconverting signals between the 1300 nm and 



980 nm wavelength bands. We note the qualitative dif- 
ference with respect to previous studies using microdisk 
optomechanical resonators, which focused on topics like 
sensitive optical transduction of motion [15j, including 
at GHz frequencies [lj|, mechanical oscillation and cool- 
ing [13], and optical mode tuning [lij . 

Microdisks are fabricated in 350 nm thick SisN4- 
on-silicon, as described in the Supplementary Mate- 
rial [19(. Similar SislNU microdisks have exhibited k/27t > 
100 MHz [2(1, so that a 10 ^ui disk diameter was chosen 
since it is small enough to place the mechanical radial 
breathing mode (RBM) frequency above the anticipated 
K, and large enough to avoid optical radiation loss due to 
imperfect total internal reflection. Finite element method 
simulations of the mechanical modes further inform our 
choice of geometry [ljl . Low mechanical dissipation re- 
quires the coupling between the breathing mode and lossy 
modes of the supporting silicon pedestal to be minimized. 
Assuming a conical frustrum for the pedestal shape, we 
calculate mechanical modes as a function of pedestal ra- 
dius at the interface with SisN4 (Fig.QJe)), with displace- 
ment amplitude profiles shown for two modes at a radius 
of 500 nm. We also plot the mechanical mode frequencies 
and Qs due to clamping losses, following the approach of 
Anctsberger et al. [2l|. We find that the pedestal radius 
should shrink to < 100 nm to achieve optimum Q; in 
practice, we have been able to achieve values of about 
100 nm, where the theoretical value is Q ps 4xl0 4 . 

We first characterize basic optical and mechanical 
properties of the system (Fig. [TJa)). Light is coupled 
to the devices using an optical fiber taper waveguide 
(FTW), and optical modes are measured by sweeping a 
1300 nm or 980 nm band tunable diode laser and record- 
ing the transmitted signal through the FTW on a 1 MHz 
bandwidth photodector. In the 1300 nm band, only 
transverse electric (TE) polarized modes are observed, 
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FIG. 1. (a) Measurement setup and scanning electron mi- 
crograph of a representative device, (b) Room-temperature 
1300 nm and 980 nm band optical modes, (c) Schematic 
of sideband spectroscopy and the microdisk. a cw /a ccw are 
amplitudes for clockwise/counterclockwise (cw/ccw) optical 
modes, which couple to the radial breathing mode (amplitude 
a;) at a rate go that is parametrically-enhanced to a rate G 
through a control field that injects N photons into the cavity. 
Ki and re e are the intrinsic loss and waveguide-cavity coupling 
rates, respectively, 7 m is the mechanical damping rate, and 
j0 is a backscattering rate that couples the cw/ccw modes. 
Pz/Pr/Pt are the incident, reflected, and transmitted pow- 
ers through the FTW. (d) Room-temperature mechanical Q 
as a function of input optical power (laser blue-detuned from 
the optical mode) for the radial breathing mode. The inset 
shows the mechanical mode spectrum at the three highlighted 
optical powers. Error bars are from the 95 % confidence inter- 
vals in Lorentzian fits of the spectra [lij.(e) Finite-element- 
method calculated frequency (top) and mechanical Q (mid- 
dle) for two modes, labeled @ and (B), as a function of top 
pedestal radius. The modes anti-cross at a radius of 200 nm. 
(bottom) Displacement profiles for a radius of 500 nm. 

while in the 980 nm band, both TE and TM (trans- 
verse magnetic) modes are seen. Mechanical modes are 
measured at low optical power (w 0.2 /iW) by tuning 
the laser wavelength to the shoulder of an optical mode, 
so that the fluctuating optical mode wavelength induced 
by thermal-noise-driven mechanical motion is converted 
to intensity modulation [![. This optical signal is de- 
tected with a 1 GHz bandwidth avalanche photodiode 
(APD) and resolved on a real-time electronic spectrum 
analyzer. At room temperature and under moderate vac- 
uum (0.002 Pa w 2xl0~ 5 torr), Q as high as 2xl0 6 at 
both 1300 nm and 980 nm are measured, on par with pre- 



vious work on SiaN4 microdisks [2(J]. Figure [ljb) shows 
a representative scan for the device we focus on, whose 
total (intrinsic) Q a is 1.5xl0 6 (1.7xl0 6 ) at 1308 nm, and 
6.7xl0 5 (7.6xl0 5 ) at 990 nm. While many mechanical 
modes are observed, we focus on the fundamental radial 
breathing mode at Lo m /2n w 625 MHz, whose mechani- 
cal Q (Q m ) is 6xl0 3 (inset to Fig. Old)). Increasing the 
optical power with the laser blue-detuned with respect 
to the cavity mode results in a pronounced narrowing 
of the mechanical mode by over an order of magnitude, 
as the laser drives the system into regenerative oscilla- 
tions [1]. On the other hand, strong pumping on the 
red-detuned side of the cavity is limited by instability 
due to thermo-optic effects. To improve this, the device 
is cooled in a liquid He flow cryostat to rj 10 K. Gener- 
ally, we see an increase in the undriven Q m , to lxlO 4 , 
and some degradation of Q (by as much as a factor of 3), 
which we attribute to cryo-gettering of material on the 
sample. The splitting in the doublet mode optical trans- 
mission spectra, produced by backscattering that couples 
the traveling wave modes of the disk into frequency- split 
standing wave modes 22j, also changes. These changes 



vary from cooldown to cooldown. 

In the context of cavity optomechanics, EIT is an in- 
crease in the transmission of a near-resonant probe beam 
(frequency U) p ) through an optical cavity (at w ) that re- 
sults when a strong control field (at lj c ) is red-detuned 
by A oc = u Q — cj c = u> m (Fig. [He)), so that anti-Stokcs 
photons generated by scattering of the control field by 
the mechanical resonator interfere with the probe and 
create a transparency window for it 0, HI HH . We focus 
on the reflected signal from the cavity, present in our mi- 
crodisks due to the aforementioned backscattering, and 
for which EIT results in a narrow dip in the reflection 
peak. The probe is derived from the control field using 
an electro-optic amplitude modulator (EOM) to generate 
higher and lower frequency sidebands, only one of which 
is coupled into the cavity since it operates in the resolved 
sideband limit. The EOM is driven by a network ana- 
lyzer, with the modulation frequency swept to vary the 
probe-control beam detuning A pc = uj p — lo c . The probe 
and control fields pass through a circulator before going 
to the device, whose reflected signal is demodulated by 
the network analyzer. A wavemeter tracks the control 
field wavelength and is used in feedback to lock A oc [lj| . 

We first study the system in the 1300 nm band, by 
sweeping A pc at different values of A oc . For each A oc , 
we take two sweeps of A pc , over a broad range to trace 
the overall cavity reflection spectrum, and over a narrow 
range to resolve the narrow EIT dip [l9( . A series of spec- 
tra is shown in Fig.[2ja), where Q has slightly degraded 
compared to its room temperature value, with k/2it s» 
180 MHz. As observed in previous work H 5], the EIT 
dip always appears at A pc =w m , and its width is given 
by the total mechanical damping rate 7, which is the 
sum of the intrinsic rate 7 TO and the optomechanically- 
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FIG. 2. EIT/EIA measurements for a (a)-(d) 1300 nm and (e)-(g) 980 nm mode, (a) Reflection spectra with the optical 
mode-control field detuning (A oc ) stepped and the probe-control detuning (A pc ) swept. Blue (red) curves are taken over a 
broad (narrow) range of A pc . (b) Zoomed-in reflection spectra (open circles are data, solid lines are fits) and (c) Depth of the 
EIT dip as a function of control field power, for A oc ~ w m . The upper inset shows the broad reflection spectrum corresponding 
to the data in (b)-(c), for which Q had degraded with respect to that in (a). The lower inset shows the width of the EIT dip 
in (b) as a function of control field power. Determination of the error bars is described in the Supplementary Information, (d) 
Zoomed-in reflection spectra for A oc = —0J m , for which EIA is observed, (e) Broad, and (f)-(g) narrow reflection spectra as a 
function of A pc for a cavity mode in the 980 nm band. In (e)-(f), A oc w cj m , while in (g), A oc ~ -w m 



induced damping rate 7om 
operativity parameter. C = 



= C'jm, where C is the co- 
|g, where G = g VN is 
the parametrically-enhanced coupling rate provided by 
N intracavity control field photons. We reach C ~ 0.45 
before the system is thermally unstable. 

To more quantitatively understand the EIT effect and 
assess the zero-point optomechanical coupling rate go , 
we set A oc w cj m and measure the normalized re- 
flection spectrum as a function of control field power 
(Fig. GJb)). For this set of measurements, Q degraded 
upon cooldown to k/2it m 550 MHz (inset to Fig. E^c)), 
while the mechanical mode had j m /2ir « 63 kHz. We 
fit each normalized reflection spectrum using the theory 
presented in Ref. [H, where the normalized amplitude re- 
flection coefficient r(A pc ) is given as: 

r ( A p c ) = : . 2«(A oc -A pc f ~ W 



1 



-+1 



i?(A pc ) = |r| 2 is the normalized reflected intensity, and in 
the limit of A oc =w m , AR = 1 - 1/(1 + C) 2 . The con- 
trast in the EIT dip, AR = 1 — i?(A pc = uj m ), is then 
plotted as a function of control field power in Fig. [2fc). 



Recalling that C 



K7„ 



, we use Eq. ([TJ in a fit to ex- 



tract go/2vr s» 7.8 kHz, where knowledge of the FTW loss, 
cavity mode transmission spectrvrm, and A oc allows us to 
estimate N at each pump power . We see good agree- 
ment between experiment and theory, further supported 
by the linear increase in EIT dip width with pump power 
(lower inset to Fig. [2jc)). However, the experimentally 
determined go is nearly two times smaller than that pre- 
dicted from finite clement simulations (<?o/27r 15 kHz) 
that exclusively consider the contribution due to moving 
dielectric boundaries. A study of whether photo-elastic 



effects may play a role [24| in this discrepancy is un- 
derway. Finally, shifting the control field to the blue- 
detuned side of the optical mode (A oc ps —u> m ) results in 
electromagnetically induced absorption (EIA), evidenced 
by a peak in the center of the reflection spectrum that 
increases with optical power (Fig. |2Jd)). 

Whispering gallery mode cavities can support high-Q 
optical modes over a broad range of wavelengths. As a 
result, multiple, broadly spaced optical modes may be ex- 
pected to couple to the same mechanical mode, which we 
observe in our measurements. Figure [UJe)-(g) presents 
EIT and EIA data for a 980 nm band optical mode 
(k/2tt ps 750 MHz) that couples to the 625 MHz radial 
breathing mode. We observe similar behavior as seen for 
the 1300 nm band mode, though the maximum contrast 
of the EIT dip and EIA peak are a little smaller, likely 
due to lower Q Q which is not fully compensated by higher 
optical power before thermal instability sets in. 

One application of multiple optical modes coupled to 
the same mechanical mode is in wavelength conversion, 
as outlined in theory 12] and recently demonstrated 

and silica micro- 



13] 



in silicon optomechanical crystals 
spheres 14J. The application of two control fields, each 
red-detuned from an optical mode by u m , opens up a pair 
of transparency windows over which wavelength conver- 
sion, mediated by the mechanical resonator, can occur. 
Input signals over a bandwidth set by the damped me- 
chanical resonator can be upconverted or downconverted, 
with an internal conversion efficiency that depends on the 
cooperativity achieved for each optical mode. 

We proceed following the recent experiments of Hill 



et al. 131] . as shown in the setup in Fig.lBTa). RF switches 



are used to alternate between EIT spectroscopy on the 
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FIG. 3. Wide-band wavelength con- 
version experiment: (a)-(b) Experi- 
mental schematics. As in Fig. [lja), 
EIT spectroscopy in the 1300 nm 
and 980 nm bands is performed on 
the reflected signal from the cav- 
ity, with fast switching between the 
two bands enabled by RF switches. 
In frequency upconversion (down- 
conversion), the 1300 nm (980 nm) 
laser is modulated to generate an 
input probe field, which is con- 
verted to the 980 nm (1300 nm) 
band through application of the 
980 nm (1300 nm) control field. 
The frequency converted field trans- 
mitted past the cavity is measured 
on the spectrum analyzer. Along 
with schematics, part (b) shows the 
optical mode transmission spectra 
for the measured device, (c) Up- 
converted (top) and downconverted 
(bottom) signals measured by the 
spectrum analyzer, measured as the 
relevant probe-control field detun- 
ing (A pc ) is swept and the 990 nm 
(top) or 1300 nm (bottom) control 
field power is stepped. 



reflected signal from the cavity in the 1300 nm band and 
in the 980 nm band. These EIT measurements are done 
to ascertain A oc for each mode, as well as the cooperativ- 
ity achieved. For the experiments that follow, A oc w u> m 
for both control fields. Focusing first on frequency up- 
conversion, we amplitude modulate the 1300 nm laser to 
generate an input probe signal that is detuned from the 
control field by A pc . The control field and input probe 
signal in the 1300 nm band are combined with the con- 
trol field in the 980 nm band using a wavelength division 
multiplexer (WDM), and sent into the FTW-coupled mi- 
crodisk. Light exiting the device is separated into the 
980 nm and 1300 nm bands using a WDM and bandpass 
filters (BPFs). Since the converted tone in the 980 nm 
band is detuned with respect to the 980 nm band con- 
trol field by A oc , the beating of the two fields is recorded 
on an electronic spectrum analyzer after detection by a 
1 GHz APD. The input probe-control field detuning A pc 
is then swept in order to assess the conversion bandwidth. 
For frequency downconversion, the roles of the 980 nm 
band and 1300 nm band lasers are reversed (Fig. EKb)). 

Figure [3Jb) shows the optical transmission spectrum 
for modes in the 1300 nm and 980 nm bands. In com- 
parison to the data from Fig. [T] and [2] the Q a here are 
lower, with k/2it w 1 GHz and 3 GHz in the 1300 nm 
and 980 nm bands, respectively. The increased losses are 
both due to additional sample degradation and the need 
to simultaneously achieve a reasonable coupling level to 
both optical modes (in Fig. [5J this was individually op- 



timized). Figure (3Jc) shows a series of converted spec- 
tra for frequency upconversion (top) and downconver- 
sion (bottom), where each trace is acquired as A pc is 
swept. For upconversion (downconversion), the 1300 nm 
(980 nm) control field power is held fixed, while the 
980 nm (1300 nm) control field power is stepped. Both 
graphs show an increase in the converted signal strength 
with control field power, though internal (external) con- 
version efficiencies [l9[ are limited to ~ 20.9 % (0.63 %). 

The maximum expected internal conversion efficiency 
is estimated from the cooperativities of the two modes, 
Ci and C 2 USHl, as rj = 4CiC 2 /(l + C x + C 2 ) 2 . While 
in Fig. [21 Ci=C 2 =0.45 was reached, corresponding to a 
maximum r} = 22%, the lower Q s measured here caused 
a reduction in C\ and C 2 that was not fully compensated 
by additional optical power, thus explaining the lower 
efficiency. The external conversion efficiency of 0.63% is 
given by the product of the internal efficiency and the 
waveguide-cavity incoupling and outcoupling efficiencies, 
which are 29.2 % and 10.3 %, respectively. 

Eventually, one might envision using such wavelength 
conversion in connecting quantum optical technology 
at wavelengths below 1000 nm (e.g., quantum memo- 
ries, Si single photon counters) and at 1300 nm (low 
loss/dispersion optical fibers). Several improvements 
must be made for such quantum frequency conversion (2o| 
to be feasible. As already described, higher coopera- 
tivities and more efficient waveguide-cavity coupling are 
needed to improve the conversion efficiency. The gen- 
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eration of noise photons is also a crucial consideration 
for applications such as frequency conversion of single 
photon states [26j], and has contributions that stem from 
the limited sideband resolution of the system (e.g., from 
Stokes scattered photons) and the thermal occupancy of 
the mechanical resonator ((n) t h ~ 330 for the mechanical 
mode at 10 K; laser cooling in our experiments reduces 
this by at most ~ 33 %). In contrast, recent work in 
optomechanical crystals [13| cooled the system close to 
the ground state {(n)th ~ 3), whereas experiments done 
at room temperature in microspheres |14| have an even 
larger noise contribution ((n)th ~ 4xl0 4 ). 

Reduced mechanical damping is one key to improved 
device performance. While Q m > 10 6 has been demon- 
strated for MHz frequency modes in S13N4 [27|, to our 
knowledge, such Qs have not yet been shown at near- 
GHz frequencies. Recently, a frequency-Q m product of 
2xl0 13 Hz was demonstrated at room temperature for 
10 MHz modes of SiaN4 membranes j28|; our correspond- 
ing value is 6xl0 12 Hz, suggesting that improved Q m 
may be possible. As simulations (Fig.QJe)) indicate that 
coupling to pedestal modes may be a limitation on Q m , 
smaller pedestal diameters or adjusting the geometry to 
isolate the periphery of the disk from a central, pedestal- 
supported region [21| may provide improvements. 

SisN4 devices may benefit by supporting a large intra- 
cavity photon number iV without nonlinear loss which, 
for example, influences Si optomechanical crystals at 
N > 300 0]. We are able to reach N w 6xl0 4 , and 
further increases are not limited by nonlinear loss, but 
instead the thermal stability of the device. Techniques to 
lock A oc beyond the simple (and relatively slow) proce- 
dure we have adopted [19( would be of significant benefit. 
Finally, we note that increased 170 could dramatically im- 
prove performance, given the squared dependence of the 
cooperativity on this parameter. While large increases in 
go are unlikely for our microdisks, since further decreases 
in diameter will begin to result in radiation losses, recent 
designs of SisN 4 slot mode optomechanical crystals [III 
suggest that such a system may be able to combine high 
frequency, large go, and large N. 
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SUPPLEMENTAL INFORMATION 
FINITE ELEMENT SIMULATIONS 

We performed finite element method simulations of 
our microdisk structures to gain a clearer physical un- 
derstanding of the optomechanical system. 

Obtaining optical and mechanical cavity modes 

Our microdisk cavities support whispering-gallery op- 
tical modes (WGMs) of the form E(r, z) exp (irruft), where 
E(r) is the electric field distribution on the rz plane (in 
cylindrical spatial coordinates) and m is the azimuthal 
order. Such WGMs were calculated by solving the 2D 
axially-symmctric electromagnetic wave equation 

V x V x E = e(r) (^) E. (2) 

via the finite element method, with a formulation that 
employed respectively edge and node elements for the 
transverse and longitudinal electric field components, 
and perfectly-matched layers to simulate open bound- 
aries. Solving the eigenvalue problem in eq. (j2J for a spe- 
cific azimuthal order m produced eigenmodes with com- 
plex frequencies Lu m and optical quality factors Q m = 
Re{w m }/(2Im{o; m }). Because WGM fields are mostly 
concentrated in the periphery of the resonator, interac- 
tion with the Si pedestal that supports the SiN microdisk 
is negligible, and thus completely ignored in the optical 
mode calculation. 

For disk radius R ~ 5 fim and SiN thickness t « 
350 nm, modes with radiation-limited quality factors ex- 
ceeding 10 8 can be found both in the 980 nm and 1300 nm 
bands. Table Q] shows wavelengths and quality factors of 
the calculated TE and TM polarized WGMs, which cor- 
respond well with. The good agreement between the ex- 
perimental and calculated values was achieved by tuning 
the parameters of the cavity within reasonable bounds 
(thickness t = 340 nm and n = 1.99 for the refractive 
index of SiN), compared to experimentally determined 
values. 

The mechanical modes of the microdisk structure were 
obtained by solving the equation of motion for the 
displacement field Q(r), assuming anisotropic materi- 
als (Io| . In our simulations, the Si pedestal supporting 
the SiN microdisk was represented by a conical frustrum 



TABLE I. Calculated and experimental whispering gallery 
modes 



Afem (nm) 


Qo,FEM X 10^ 


Polarization 


Aexp. ( nm ) 


971.28 


10 


TE 


970 


987.37 


4.5 


TE 


984 


975.36 


5.7 


TM 


974 


989.61 


2.4 


TE 


990 


1281.23 


1.0 


TE 


1280 


1309.04 


0.5 


TE 


1308 



whose angle and height were equivalent to those deter- 
mined from scanning electron microscope images of fab- 
ricated structures. A zero displacement boundary con- 
dition was enforced at the bottom surface of the conical 
frustrum, corresponding to the region where the pedestal 
meets the Si substrate. 

For the 10 fim diameter disks considered, a first-order 
radial breathing mode (RBM) is obtained in the vicinity 
of / ro =625 MHz. The displacement profile of the RBM, 
as shown at the bottom of Fig. 1(e) (mode (§)), is primar- 
ily in the radial direction and confined to the disk plane, 
with relatively small vertical displacement in the disk- 
pedestal contact area. Due to its azimuthal symmetry, 
the RBM is expected to display preferential optomechan- 
ical coupling to the microdisk's WGMs. 

Optomechanical coupling 

The shift in the frequency f a of a particular opti- 
cal resonance due to displacement of the nanostructure 
boundaries produced by a mechanical resonance at fre- 
quency f m is quantified by the optomechanical coupling 
9om = doj /dx = w /ioM; here, x is the cavity boundary 
displacement, tu = 2nf and Lom is an effective optome- 
chanical interaction length [l(| ■ The effective length Lom 
can be estimated via the perturbative expression 

r 2/dy £ |E| 2 

JdA(Q-n) (Ae|E||| 2 -A(e-i)|DJ 2 ) 

Here, E and D are the modal electric and electric dis- 
placement fields, respectively, Ae = €diei.—£air, A(e _1 ) = 
e diei. ~ e ain an d ^diel.,air are the permittivities of the mi- 
crodisk material and air. The mass displacement due to 
the mechanical resonance is given by Q, and the nor- 
mal surface displacement at the structure boundaries is 
Q • n, where n is the surface normal. The integral in the 
denominator is performed over the entire surface of the 
nanostructure. 

The optomechanical coupling g om can be converted 
into a pure coupling rate g between the optical and me- 
chanical resonances: go — x zp f ■ g m, where x zp f — 
y/h/2muj m is the zero point fluctuation amplitude for 
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TABLE II. Optomechanical coupling parameters between 980 
nm and 1300 nm TE-polarized WGMs and 625 MHz RBM 



A (nm) 


Qo x 10 8 


g om /27r (GHz/nm) 


g /2ix (kHz) 




988.41 


4.8 


42.7 


19.9 


7.10 


1281.80 


1.0 


33.2 


15.5 


7.05 



mechanical displacement and m is the motional mass 
of the mechanical resonance at frequency u> m . The 
motional mass can be obtained from the displacement 
Q and the nanobeam material density p by m = 

For the RBM at f m rj 625 MHz above, m rj 60 pg 
and x zp f ss 0.5 fm. Table [XT] shows calculated values of 
L om , 9om an d the vacuum coupling rate go between the 
RBM and TE-polarized WGMs in both the 980 nm and 
1300 nm bands. The calculated g om values are smaller 
than those estimated via the expression goM — w /R, 
with w the optical mode frequency and R the disk ra- 
dius, which is commonly employed for the case of large 
disks. For smaller radius disks such as considered here, 
optical fields extend considerably into the vacuum re- 
gions surrounding the dielectric, resulting in an effective 
optomechanical length L om larger than the actual phys- 
ical dimension R, as seen in Table \H\ 



Clamping losses 

To estimate mechanical clamping losses of the me- 
chanical resonances, we adopted the method of ref. 21| . 
in which the contact region between the disk and the 
pedestal is modeled as a membrane that radiates acous- 
tic energy with power 



cL4|Q(r).z| 2 . 



(4) 



where p is the density of the disk material and c is the 
speed of sound, fi m is the angular frequency of the me- 
chanical resonance and |Q(r) • z| 2 is the out-of-disk-plane 
displacement over the contact area A. If W mec h is the 
stored mechanical energy of the resonance, the latter's 
mechanical quality factor can be estimated as 



P 



n m w„ 



(5) 



It is apparent in Fig.l (e), that, for large top pedestal 
radius, the mechanical quality factor for the RBM tends 
to decrease with increasing pedestal radius, a result of 
the increased contact area through which energy may 
radiate. At the same time, for a contact areas with 
radius near 200 nm, the RBM mixes with a secondary 
('pedestal') mode that displays a large vertical displace- 
ment at the disk center, as observed in the profile of mode 



@. Such mixing is evidenced by an anti-crossing between 
the green and blue mechanical frequency curves in Fig. 
1(e), and by a steep decrease in quality factor in the 
neighborhood of the anti-crossing. Reduced quality fac- 
tors in this range are associated with the relatively large 
vertical displacement of the pedestal mode. For radii 
below 200 nm the RBM quality factor increases again 
(green curve in Fig. 1(e)), as the pedestal mode is driven 
towards lower frequencies. 

While the geometrical details of the Si pedestal must 
be included for an exact determination of the pedestal 
mode frequencies, we have observed that, generally, 
smaller pedestal top radii were necessary to ensure that 
such modes would be out of the range of the RBM. Our 
highest mechanical quality factors were measured for a 
pedestal radius of 100 nm, which is essentially the small- 
est size that we can produce with our current fabrication 
process. 



DEVICE FABRICATION 

The device fabrication started with a bare silicon 
wafer, as shown in Fig. [4] A layer of 350 nm thick stoi- 
chiometric SisN4 film was grown by low pressure chem- 
ical vapor deposition (LPCVD), with a process- induced 
internal tensile stress of ~ 800 MPa, as measured by 
the wafer bowing method. A 500 nm thick positive- 
tone electron beam (E-beam) resist was spin-coated on 
the SisN4 film, followed by E-beam lithography and de- 
velopment in hexyl acetate at 8 °C. The patterns were 
then transferred into the SiaN4 layer by an 02/CHF 3 /Ar 
inductively-coupled plasma reactive ion etch (RIE) . This 
RIE step apparently leaves a thin layer of Si02 on top of 
the exposed Si surface, preventing any subsequent KOH 
undercut (Fig. 0|b)). To rectify this, before the e-beam 
resist was removed, an additional SF6/C4F8 inductively- 
coupled plasma RIE was carried out to remove the Si02 
layer (as well as some Si), while the residual E-beam re- 
sist protected the S13N4 device layer. After resist removal 
using a stabilized H2SO4/H2O2 solution, the sample was 
undercut in a 20 % KOH bath, as this concentration has 
an etch rate that is relatively insensitive to temperature 
fluctuations. The KOH etch was performed in two steps 
to achieve a small (< 200 nm) top pedestal diameter un- 
der the microdisk. We started to etch the sample at 40 °C 
to quickly remove the bulk Si, with periodic inspection of 
the pedestal size under an optical microscope every 6 to 
10 minutes. When the pedestal size was less than 1 //m, 
the undercut was continued at a temperature of 20 °C to 
reduce the etch rate. The undercut was stopped every 4 
minutes for inspection, until one or more microdisks were 
completely released (etch rate variations across the chip 
prevent all devices from being completely released for the 
same undercut time). The pedestal diameters for the re- 
maining (intact) microdisks were typically <200 nm, as 
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(a) 



LPCVD Si N growth 



(b) 



E-beam lithography 
&CHF 3 /0 2 /Ardry etch 



(c) 



SF7CH dry etch 



Si 3 N 4 



Si 



. E-beam resist 
SiO„ 



stable over several scans in a row. Typical locking times 
are on the order of a few seconds (for lower powers) up to 
a couple tens of seconds (for the highest powers) . Power 
levels beyond that shown in Fig. 2 result in thermally 
unstable behavior, where the cavity mode shifts too far 
away from the control laser for the above approach to be 
effectively applied. 

Once the locking procedure kept A oc within a prede- 
termined acceptable range around the target value (e.g., 
50 MHz), EIT measurements consist of sweeps of A pc 
recorded over two frequency ranges - the first over a broad 
range to provide an overall picture of the cavity reflection 
spectrum, and the second over a narrow range to record 
the details of the EIT dip. For each set of parameters 
(optical powers and A oc ), these sweeps were repeated 16 
times and the data plotted in Fig. 2 is an average of these 
scans. 



(d) 



KOH undercut at 40 °C 



FITTING THE DATA 



(e) 



KOH undercut at 20 °C 
□ 



FIG. 4. Processing steps for fabrication of the SisN4 mi- 
crodisks. 



determined by imaging with a scanning electron micro- 
scope. 



Experimentally-measured optical transmission spec- 
tra are fit using coupled mode theory for a resonator- 
waveguide system. For single dips, this yields a 
Lorentzian function, while doublet modes are fit using 
a model that takes into account backscattering due to 
surface roughness, as described elsewhere 3l|, |32j. In the 
limit of large doublet spitting with respect to the cavity 
linewidths, a pair of Lorentzians can accurately fit the 
data. Mechanical spectra were fit with single Lorentzian 
functions. 

Fitting of the EIT signals (Fig. 2(b)) is based on Eq. 
(1) from the main text, which we repeat here: 



LOCKING OF A oc 

In EIT experiments, the microdisk cavity is heated 
when the control laser reaches high power levels, resulting 
in a red-shift of the resonances due to the thermo-optic 
effect. It is thus difficult to maintain a fixed detuning 
(A oc ) between the control field and optical cavity with- 
out feedback. To reduce this instability, we lock A oc by 
adjusting the control laser wavelength to track the red- 
shifted cavity resonance, following a procedure similar to 
that outlined in [3(J. First, a small portion of the control 
field is split off and fed to a wavemeter (Fig. 1) equipped 
with a PID controller whose output is fed back to the 
laser to fix its absolute wavelength at a specified value. 
This specified value is determined by sweeping A pc over 
a broad frequency range (0 to 3 GHz) and fitting the ac- 
quired data by the vector network analyzer to determine 
the real value of A oc . The control field wavelength is then 
adjusted to ensure that the desired A oc is achieved, and 
the process repeated until the fit value of A oc remains 



r(A pc 



2»(A oc -A pc ) 



(6) 



2i(^ m -A pc 
7m 



The values of k, A pc , and A oc were acquired from the 
broad scans of the EIT signal. The intracavity photon 
number is determined as: 



N = -^—y^ATQi( — 

nuj u> 



(7) 



where h is the Planck constant divided by 27r, £ is the 
FTW transmission, AT is the depth of the optical res- 
onance in the transmission spectrum (AT=1 at critical 
coupling), Qi the intrinsic optical Q, and Pi n is the opti- 
cal power at the FTW input. We measured a background 
signal with the FTW far away from the device to account 
for the frequency-dependent response of the APD, EOM, 
and network analyzer. The fitting was carried out on the 
experimentally acquired EIT data divided by the back- 
ground signal. 
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626 626.5 627 627.5 626 626.5 627 627.5 

frequency (MHz) < b * frequency (MHz) 



FIG. 5. (a) Raw data (blue) and the fitted curve (red) of 
the 1300 nm EIT signal at the control power of 21.9 fiW. (b) 
Difference between the raw data and the fitted curve shown 
in (a). 

The fitting of the EIT width versus control power (inset 
in Fig. 2(c)) is based on the expression for the total 
mechanical damping rate 7 = 7 m (l + C). The fitted 
value of go from the EIT width is approximate equal to 
that obtained from the EIT depth fitting, with an error 
of sa 6 %. 



UNCERTAINTY ESTIMATES 



The error bars in Fig. 2(c) arise from the noise of 
recorded raw EIT data. By subtracting the fitted curve 
from the raw EIT data (Fig. Ufa)), we obtained the noise 
of the EIT signal (Fig.JSJb)). The one standard deviation 
value of this noise signal is considered as the uncertainty 
of the EIT signal depth and is plotted as the error bar in 
Fig. 2(c). By combining the uncertainty of the reflected 
EIT signal and the slopes of the fitted curves at both the 
left and right shoulders, we obtained the uncertainty of 
the width of the EIT signal. This width uncertainty is 
plotted as the error bars in the inset figure of EIT width 
versus the optical power in Fig. 2(c). 

The uncertainty in Q m shown in Fig. 1(d) is given by 
the 95 % confidence intervals of the Lorentzian curve fit- 
ting of the experimentally measured mechanical spectra. 



